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Abstract: We present a model for the formation of the charmed pentaquark @ c in hard 
scattering processes such as deep inelastic scattering, e + e~ annihilation, and high-energy 
pp collisions. The model assumes that the cross section for C formation is proportional 
to the rate of production of pD*~ (or pD* + ) pairs in close proximity both in momentum 
space and in coordinate space. The constant of proportionality is determined from the G c 
cross section in deep inelastic scattering as reported by the HI experiment. The HERWIG 
Monte Carlo is used to generate simulated DIS events and also to model the space-time 
structure of the final state. Requiring the proton and the D* be within a 100 MeV mass 
window and separated by a spacelike distance of no more than 2 fm, we find that a large 
"coalescence enhancement factor" F co ~ 10 is required to account for the HI signal. The 
same approach is then applied in order to estimate the number and characteristics of C 
events produced at LEP and the Tevatron. 
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1. Introduction 

During the last year, many experiments have reported the observation of an exotic uudds 
pentaquark baryon 0+ [1, 2]. Shortly after the initial experimental reports, Ref. [3] pos- 
tulated the existence of an anti-charmed analogue C with quark content uuddc, and com- 
puted its likely properties, using the diquark-triquark configuration earlier conjectured to 
describe + [4]. 

Ref. [3] suggested that C is an isosinglet with J p = i + and estimated its mass 
at 2985 ± 50 MeV. It also discussed another possible exotic baryon resonance containing 
heavy quarks, the 6^, a uuddb state, and estimated M(B^) = 6398 ± 50 MeV, noting 
that these states should appear as unexpectedly narrow peaks in D~p, D°n, B°p and 
B + n mass distributions. Possible search opportunities in e + e~ have also been suggested 
in Refs. [5, 6, 7]. 

In March, HI reported [8] observing a narrow resonance in D*~p and D* + p invariant 
mass combinations in inelastic ep collisions at Ecu of 300 GeV and 320 GeV at HERA. 
The resonance has a mass of 3099 ±3 (stat.) ± 5 (syst.) MeV and a measured Gaussian 
width of 12 ±3 (stat.) MeV, compatible with the experimental resolution. The experiment 
interpreted this resonance as an anti-charmed baryon with a minimal constituent quark 
composition of uuddc, together with the charge conjugate. 

A parallel analysis by ZEUS sees no signal [9]. Null results have also been reported 
in papers and unpublished reports at conferences by FOCUS [10], ALEPH [11] and CDF 



Since HI observes the resonance in both D*~p nd D* + p, it is extremely unlikely that 
their result is due to a statistical fluctuation. 

Since the ZEUS and HI experiments are quite similar, a possible interpretation of 
these results is that O c does not exist and the HI observation is due to some yet to 



[12]. 
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be understood systematic effect. The discrepancy between HI and ZEUS is of a purely 
experimental nature. It is now being carefully examined by the two collaborations. 

Putting the ZEUS vs HI issue aside for the time being, an important question is: 
what are the implications of the null results in non-HERA experiments? Since production 
mechanisms might be different in different experiments, an a priori possible interpretation 
is that these experiments have not accumulated enough data or that more sophisticated cuts 
need to be applied in order to separate the likely small C signal from the background. In 
the present paper we examine the quantitative implications of a possible Hl-like production 
mechanism at LEP and the Tevatron. 

One can try to apply a simple counting argument to the HI data in order to estimate the 
expected number of events at LEP. The HI data sample of 75 pb _1 of DIS events contains 
roughly 1.5 x 10 6 events with charm. This should be compared with the corresponding 
numbers for Z° decays at LEP: 

r(hadrons)/r total = 0.70, 
r(cc)/r(hadrons) = 0.17, (1.1) 

so T(cc)/T^ ^ a Y = 0.12. After cuts, each of the LEP experiments has several million 
hadronic Z° decays. Specifically, ALEPH [11] has 3.5 xlO 6 , corresponding to approximately 
6 x 10 5 charm events. 

We can use this number to make a very crude estimate of the number of expected 
anticharm pentaquark events at LEP. If the production mechanisms at HI and LEP are 
roughly similar, the expected number of anticharmed pentaquarks events at each of the 
four LEP experiments is about 40% of the number of such events at HI. 

In the rest of the paper we carry out this estimate in more detail, with a coalescence 
model providing a specific possible mechanism for 6 C formation from D*p, and HERWIG 
generating the corresponding D* and p distributions. 

The result of this analysis is that each of the LEP experiments should have observed 
between 25 and 40 charmed pentaquark events of the HI type. This takes into account 
typical experimental cuts and efficiencies. 

At the Tevatron, the integrated luminosity already collected in Run II is of the order 
of 200 pb _1 . If the production mechanism is the same as in HI, we conclude that even 
with a very low overall detection efficiency of only a few per mille, this would provide an 
unmissable signal on the order of 10 5 events. 

In the following sections of the paper we develop the model, extract the relevant coa- 
lescence parameter from the HI data and then apply it to LEP and Tevatron experiments. 
We then discuss several possible interpretations of these results. 

2. Coalescence model 

We construct a simple but, we believe, physically appealing model to attempt a quantitative 
estimate of the expected production rate of O c via coalescence of a D* and a nucleon. 
We use the Monte Carlo program HERWIG [13, 14] to generate simulated DIS events with 
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HERA kinematics, selecting those which include either p and D*~ or p and D* + . In addition 
to the momentum-space structure of final states, HERWIG generates an approximate space- 
time structure, which enables us to model the joint distribution of invariant mass and space- 
time separation of the pD* pairs. Assuming resonance formation through coalescence, we 
suppose that the cross section for 6 C production is proportional to the number of pD* pairs 
that have invariant masses in a small interval Am around the resonance mass and are at 
the same time formed within a small spacelike separation Ax. By adjusting the constant 
of proportionality to match the Q c signal observed in DIS by HI, we predict the G c signal 
that should be seen in other processes with a similar production mechanism. 

Before going into a detailed discussion of the model's implementation, a word of warn- 
ing is in order. The present work presents a new type of approach, going beyond the 
standard coalescence model. This is made possible by the fact that HERWIG generates 
a space-time distribution of hadrons produced in high-energy collisions. As with any new 
approach, one would like to first test it against other similar cases where the answer is 
known. This can be done to some extent by looking at (anti)-deuteron production at 
HERA and at LEP, as discussed at the end of Sec. 2.3. But the experimental situation 
with deuterons is puzzling on its own, and the analogy can be only pushed so far. As 
for other well-known resonances, their dominant production mechanism is hadronization 
of quark and gluon jets, rather than through coalescence. The relevant yields of such res- 
onances are built into HERWIG from the beginning and closely follow the experimental 
data. However, the charmed pentaquark, and especially the baffling experimental situa- 
tion, is sufficiently interesting that we believe it provides strong motivation for the study 
of the production mechanism, even though the model cannot be directly tested on other 
well-known resonances. 

2.1 Model for space-time structure 

Since the space-time structure of the final state is a crucial ingredient of the proposed 
coalescence mechanism, we start with a brief summary of the HERWIG model for this 
structure. Fuller details may be found in ref. [13]. 

During development of the final state, HERWIG assigns to each unstable object with 
four-momentum q an exponentially distributed proper lifetime, with mean value r given 



where M and T are the nominal mass and width. This formula interpolates between 
r = h/T for an object that is on mass-shell and r = fr\/~q^/(q 2 — M 2 ) for one that is 
far off mass-shell. Following sequential electroweak decays and QCD parton showers, the 
resulting space-time separations are combined to give the locations of partons at the start 
of the hadronization process. 

For light quarks and gluons, whose natural widths are small, the above prescription 
can lead to unreasonably large distances being generated in the final, low-virtuality steps 
of parton showering. To avoid this they are given a width V = VMIN2/M, the HERWIG 
parameter VMIN2 (default value 0.1 GeV 2 ) acting as a lower limit on a parton's virtuality. 



by [15] 



r = 




(2.1) 



^{q 2 - M 2 ) 2 + {Tq 2 /M) 2 
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Hadronization in HERWIG takes place through non-perturbative gluon splitting, g — > 
qq, followed by colour-singlet qq cluster formation and decay. The production vertex of 
a cluster is taken as the midpoint of a line perpendicular to the cluster's direction of 
travel, with its two ends on the trajectories of the constituent quark and antiquark. The 
production positions of primary hadrons from cluster decays are smeared, relative to the 
cluster position, according to a Gaussian distribution of width h/M where M is the mass 
of the cluster. 



2.2 DIS results 



We first use HERWIG to simulate DIS in e + p collisions with beam momenta 27.6 x 920 
GeV/c and integrated luminosity 75 pb _1 , which approximates the data sample used in 
the HI analysis [8]. For greater efficiency, we limit the target partons to charmed quarks 
and antiquarks only, by setting the process code IPRDC = 9004. Using the HERWIG default 
(MRST LO [16]) parton distribution functions, we find a total leading-order DIS-on-charm 
cross section of 18.4 nb for momentum transfers in the region 6 < Q 2 < 100 GeV 2 and 
inelasticities 0.05 < y < 0.7. We therefore generate 1.4 x 10 6 simulated events so as to 
correspond roughly to the 75 pb" 1 integrated luminosity of the HI data sample. In HERWIG 
runs with IPR0C = 9000 (DIS on all flavours of partons), we find that the contribution to 
pD* production in the region of interest from DIS on non-charmed partons is negligible. 

We use a momentum transfer range of 6 < Q 2 < 100 GeV 2 , rather than the full 
range 1 < Q 2 < 100 GeV 2 covered by the HI analysis, because the charm PDF and 
HERWIG simulation are not reliable at lower Q 2 . The predicted dependence of the D*^ 
cross section on Q 2 is shown in fig. 1 together with the recent HI data [17]. The predicted 
D* pseudorapidity distribution is shown with the HI data in fig. 2. These results can also 
be compared with ZEUS results [18] for charm production in DIS. 

The histograms in figs. 1 and 2 are for the cuts in ref. [8], 
while the data shown are for slightly different cuts (see ref. [17] 
for details), which we estimate may reduce the cross section 
by 30-40% without changing the shape significantly. In view 
of the more limited range of Q 2 accessible in the HERWIG 
simulation, a detailed comparison with the data is not possi- 
ble at present, but we judge the semi-quantitative agreement 
to be satisfactory. 

Next we select events in which both a proton and a D*~ , 
or an antiproton and a D* + , are produced. Imposing the cuts in Table 1, which approximate 
those applied in the HI analysis, but not yet selecting particular D* decays, we obtain the 
pD* mass spectrum shown by the solid histogram in fig. 3. 

The predicted joint distribution in mass and space-time separation is shown in fig. 4. 
Here we depict spacelike separations as positive and timelike as negative. There is a 
preference for spacelike separations, peaking at small values, corresponding to the formation 
of primary hadrons on a hyperboloid at a roughly universal proper hadronization time. This 
should be a general characteristic of physically reasonable hadronization models. 
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Table 1: Selection criteria 
for pD* production. 
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Figure 1: HERWIG prediction for D*^ production in DIS at Q 2 > 6 GeV 2 , together with recent 
HI data [17]. 
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Figure 2: Predicted Z?** pseudorapidity distribution in DIS, together with recent HI data [17]. 



We now suppose that the C resonance, not explicitly present in the HERWIG hadroniza- 
tion scheme, can be modelled as an enhancement in the cross section proportional to the 
number of pD* pairs with invariant masses in a window Am close to resonance and at 
the same time within the range Ax of the binding mechanism. The measures of closeness 
should be on the hadronic scale of Am ~ 100 MeV, Ax ~ 1 fm. To be specific, we choose 
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Figure 3: Predicted non-resonant pD* mass distribution in DIS. Solid: all pD* and pD* + pairs. 
Dashed: pairs with separation < Ax < 2 fm. 
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the mass interval 3050 < m < 3150 MeV and spacelike separations < Ax < 2 fm. 

The predicted mass distribution for pD* pairs with separations < Ax < 2 fm is 
shown by the dashed histogram in fig. 3. Of these pairs, 245 fall within the coalescence 
region 3050 < m < 3150 MeV. With a coalescence enhancement factor F co , we therefore 
expect of the order of 250F CO events with G c production. 

In the HI experiment, the D*^ is detected through 
the decay sequence D*^ 1 — > D ^, D° — > K~n + , 
which has a branching fraction of about 3%. Taking 
into account the HI cuts needed to identify the de- 
cay products, listed in Table 2, we estimate that this 
fraction is reduced to roughly 2%. Therefore the ex- 
pected HI G c signal is about 5F co events. Given the 
observed signal of 50 events, this implies a rather 
large coalescence enhancement factor of F co ~ 10. 

To obtain the observed pD* mass distribution from the one shown by the solid his- 
togram in fig. 3, one should multiply by the D*^ detection efficiency. The shape of the 
distribution is similar to that shown in ref. [8], but the above efficiency estimate of 2% im- 
plies an observed number of about 10 pairs per 50 MeV bin, which is considerably smaller 
than the number of about 10 pairs per 10 MeV bin found by HI. This could be due to 
an underestimate of the proton yield by HERWIG, a higher relative proton yield in charm 
production with 1 < Q 2 < 6 GeV 2 , or a large contamination of the experimentally reported 
distribution by non-protons. 

We expect the predicted proton yield in the coalescence region to be more reliable, 



Table 2: 

D* -» D° 7T S 



Selection criteria for 
D° -> K-TT+ decay. 
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Figure 4: Predicted p£>* joint mass-separation distribution in DIS. 



since additional soft contributions will tend to be more widely distributed in space-time. 
We therefore use the value F co ~ 10 throughout the rest of this paper. However, it should 
be borne in mind that an enhancement of the proton density in the coalescence region 
in DIS would imply a reduced value of the coalescence factor F co and hence a smaller 
prediction of C production in other processes. On the basis of the information available 
to us, we estimate the maximal range of values indicated by the HI signal to be 2^F co ^10. 

Of course, changing the size of the coalescence region by varying Am and/or Ax can 
be compensated by a change in F co to give the same signal. For example, one could increase 
the mass interval to 200 MeV, i.e. 3000 < m < 3200 MeV, and the maximum separation 
to 4 fm. Then a value of F co ~ 2.5 would account for the HI signal, without changing 
substantially the predictions given below. This stability against variation of Am and/or 
Ax provides a nontrivial test of the coalescence model's robustness. 

Having fixed the coalescence factor F co , we can make predictions about other aspects 
of O c production in DIS. For example, we show in fig. 5 the predicted O c production cross 
section in DIS as a function of the momentum transfer Q 2 . The rather weak dependence 
at low Q 2 helps to justify our exclusion of the region Q 2 < 6 GeV 2 . Figures 6 and 7 show 
the predicted pseudorapidity and transverse momentum distributions of the O c , computed 
from the corresponding quantities for pD* pairs in the coalescence region. 
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Figure 5: Predicted 9 C production cross section as a function of momentum transfer in DIS. 
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Figure 6: Predicted 6 C pseudorapidity distribution in DIS. 



2.3 e+e results 

The coalescence model can be used to predict the rates and distributions of 6 C production 
in other hard processes that can be generated by HERWIG, such as e + e~ annihilation. We 
simply assume that the C distributions are always in the fixed ratio F co = 10 to those of 
pD* pairs in the coalescence region of Am and Ax. 

We concentrate on the resonant reaction e + e~ — > Z° — > hadrons since this is the place 
where the LEP experiments accumulated most of the data that might show evidence of 
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Figure 7: Predicted 9 C transverse momentum distribution in DIS. 



O c production. We generate 7 x 10 5 simulated events of the process e + e~ — > Z° — > cc at 
the resonance peak, which corresponds to the approximate integrated luminosity of 140 
pb _1 collected by each experiment at or near Z° resonance. Although the process e + e" — > 
Z° — > bb is a copious indirect source of charm, the possible production of C in 6-decay 
is a separate issue and we assume that this source could be distinguished experimentally 
from direct production. As in the case of DIS, we find that, within the framework of the 
HERWIG model, the contribution to pD* production in the coalescence region from light 
quark production (e + e _ — > Z° — > qq where q = u,d, s) is negligible. 

We apply no selection criteria like those in Tables 1 and 2, since these are mostly 
irrelevant to the LEP environment. Therefore our predictions represent upper limits on 
the observable signal at LEP1 for consistency with the HI data and the coalescence model. 

Figure 8 shows the resulting pD* joint distribution in mass and space-time separation, 
to be compared with fig. 4 for DIS. The two distributions are similar, with some additional 
narrowing of the mass distribution at small separations in the case of Z° decay. 

The overall pD* mass distribution is shown by the solid histogram in fig. 9, corre- 
sponding to fig. 3 for DIS. In order to make a comparison with experiment, we need to 
take into account data reduction due to cuts and efficiencies. We generated 7 x 10 5 charm 
events, corresponding to 7 x 10 5 /0.17 = 4.1 x 10 6 hadronic events, to be compared with 
3.5 x 10 6 hadronic events reported by ALEPH [11] after cuts. ALEPH reports 1.3M proton 
candidates, whereas the proton yield is about 1 per event, so their mean proton efficiency 
is ~ 1/3. For the D* efficiency, through D*^ — > D ^, D° — > K~it + , we assume a range 
2-3%. These cut and efficiencies yield a reduction factor between 5.7 x 10~ 3 and 8.5 x 10 -3 . 
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Figure 8: Predicted pD* joint mass-separation distribution in Z a decay. 



In [11] ALEPH reports 82 D*~p events in the mass range 2.95 GeV < M(D*~p) < 
4 GeV. The total number of D~*p events in this mass range in fig. 9 is about 15,000, so 
we would expect between 85 and 130 events to be seen, in good agreement with [11]. 

The mass distribution after a cut on separation is shown by the dashed histogram in 
fig. 9. The number of pairs in the coalescence region < Ax < 2 fm, 3050 < m < 3150 
MeV is 576, roughly as expected from the 245 found in DIS with about half the integrated 
luminosity. The corresponding C production cross section on resonance, assuming a coa- 
lescence enhancement factor of F co ~ 10, is 40 pb. Taking into account the data reduction 
factors above, that implies a signal of about 25-40 events per experiment in the observed 
decay channel. 

The predicted distribution of C production with respect to the momentum fraction 
x p = 2\p\/y/s is shown in fig. 10. This distribution is similar in shape to that of the D*^, 
but slightly narrower (fig. 11). The predicted C production is strongly collimated along 
the jet axes of the final state, as shown by the pseudorapidity and transverse momentum 
distributions relative to the thrust axis (figs. 12, 13). 
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Figure 9: Predicted non-resonant pD* mass distribution in Z° decay. Solid: all pD*~ and pD* + 
pairs. Dashed: pairs with separation < Ax < 2 fm. 
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Figure 10: Predicted cross section e + e — > Z° — > 9 C ^ as a function of momentum fraction. 
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Figure 13: Predicted cross section e + e~ — > Z Q — > C X as a function of transverse momentum 
relative to the thrust axis. 

In contrast, ALEPH reports [11] no deviation from background in the relevant mass 
region, i.e. no signal. Our results then have several possible interpretations: 

(a) there is no resonance and the HI signal is due to some yet to be understood systematic 
effect; 

(b) HI observed a genuine resonance, but for some reason the production mechanism in 
DIS is different than in e + e~; 

(c) the production mechanism is the same in DIS and e + e~ but it is not described by 
the coalescence model; 

(d) some crucial ingredient in the e + e~ data analysis is not well understood. 

In this context we would like to draw the reader's attention to two examples of non- 
exotic baryon states which are definitely known to exist, but whose study seems to be 
problematic at LEP. 

The first example is A c (2880) + which is above the ND threshold. According to RPP 
[19], the only observed decay modes of A c (2880) + are A+7r + 7r~ and £ c (2550)7r. So is 
unclear why no one has seen A c (2880) + — ► DN (D*N is above threshold). 

An even more striking example is the non-observation of (anti)-deuterons at LEP. As 
already noted elsewhere [20], it is also useful to compare the rate of antideuteron and 
antiproton production in a given experiment. Such an analysis has been carried out by 
HI [21], yielding an antideuteron/antiproton ratio of d/p = 5.0 ± 1.0 ± 0.5 x 10~ 4 in 
photopro duct ion . 
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On the other hand, although the LEP experiments produced roughly one proton per Z 
decay [22] and have accumulated millions of Z° decays on tape, very little is known about 
antideuteron production at LEP. The one theoretical prediction we are aware of is Ref . [23] , 
which uses the Lund string fragmentation model to predict 5 x 10 -5 deuterons per Z° decay. 
The only relevant experimental publication we are aware of is from OPAL [24], which 
reports exactly one antideuteron candidate event which was eventually dismissed because 
it did not pass through the primary vertex. From this OPAL infers at 90% confidence 
level an upper limit on antideuteron production of 0.8 x 10 -5 antideuterons per Z° in the 
momentum range 0.35 < p < 1.1 GeV. 

A recent estimate [25] based on this data concludes that d/p < 1.6 x 10~ 4 , which is sig- 
nificantly less than the ratio reported by HI [21]. 1 The reason for this presumed difference is 
unknown at present. In analogy with O c formation, one can construct a coalescence model 
for deuteron formation. Using a similar coalescence window (m — m p — m n < 100 MeV, 
< Ax < 2 fm), the LEP upper limit implies a coalescence factor for p + n — > d, of 
Fc£* < 0.1. It would be very valuable to have more information on antideuterons from the 
LEP experiments. 

2.4 Hadroproduction results 

The © c coalescence model can also be used to predict its production in hadron-hadron 
collisions. The only limitation, as in DIS, is that HERWIG is able to model only hard 
processes. 2 In heavy quark hadroproduction the quark mass provides a hard scale, but in 
the case of charm this is not large enough to ensure perturbative reliability. We therefore set 
a minimum transverse momentum for the hard subprocess, setting the HERWIG parameter 
PTMIN to 5 GeV. This means that charmed hadron production can only be predicted at 
transverse momenta greater than this value. 

The prediction for D*^ production above 5 GeV in Tevatron Run II (pp at y/S = 2 
TeV), shown in fig. 14, is indeed in reasonable agreement with CDF preliminary data [26]. 

The model study presented here is based on a sample of 10 5 simulated charm produc- 
tion events. The HERWIG estimate of the charm cross section with the above transverse 
momentum cut is 0.15 mb. Therefore the equivalent integrated luminosity of our sample 
is quite small, around 0.7 nb -1 . 

Figure 15 shows the corresponding pD* joint distribution in mass and space-time 
separation, to be compared with figs. 4 and 8 for DIS and Z° decay, respectively. We 
see that the distribution in hadroproduction is substantially broader in both variables. 
This is due to the more complex structure of hadron-hadron collisions. In particular, the 
underlying event (due to interactions between the beam hadron remnants) often contains 
protons and antiprotons formed over a longer timescale than that required for coalescence. 



x We thank T. Sloan for discussion of this point. 

2 The results in this section include correction of an error in HERWIG in the programming of the space- 
time structure of the underlying event in hadron-hadron collisions. The quantitative changes are small. 
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Figure 14: HERWIG prediction for D*^ transverse momentum distribution in pp collisions at 
VS = 2 TeV. 




Figure 15: Predicted pD* joint mass-separation distribution in charm hadroproduction at \/S = 2 
TeV. 



The pD* mass distributions before and after a cut on separation are shown in fig. 16, 
corresponding to figs. 3 and 9. The predicted C transverse momentum distribution is 
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Figure 16: Predicted non-resonant pD* mass distribution in charm hadroproduction at \/S = 2 
TcV. Solid: all pD*~ and pD* + pairs. Dashed: pairs with separation < Ax < 2 fm. 



50 



40 



> 
o 



ft 

T3 



b 



30 — 



20 — 



10 — 



0.0 2.5 5.0 7.5 

Pt (GeV) 



10.0 



12.5 



15.0 



Figure 17: Predicted 6 C transverse momentum distribution in pp collisions at y/S = 2 TeV. 



shown in fig. 17. The number of pairs in the coalescence region < Ax < 2 fm, 3050 < 
m < 3150 MeV is 19. Taking into account coalescence factor F co ~ 10, this corresponds 
to 190 C events, i.e. to a O c production cross section of 285 nb. Since the integrated 
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luminosity already collected in Run II is of the order of 200 pb , this corresponds to 
production of some 57 million c 's. Even with an overall detection efficiency of only a few 
per mille, this would seem to provide an unmissable signal. 

3. Summary and conclusions 

We constructed a simple coalescence model for the production of the charmed pentaquark 
reported by the HI experiment. This model can also serve as a template for estimating 
resonance and bound state production cross sections in a wide range of hard processes. 

In the specific case discussed here pentaquark formation is assumed to occur when 
D*p are in close proximity in both momentum and coordinate space. Comparing the 
HERWIG-generated D*p distribution with the signal reported by HI, we determine the 
coalescence-enhancement factor F co ~ 10. This is then applied to estimate the number of 
C events in e + e~ annihilation at LEP and pp collisions at the Tevatron. For each of the 
four LEP experiments the model then predicts between 25 and 40 Hl-like C events. For 
the Tevatron a signal of many thousands of events would have been expected. 

Since both LEP and Tevatron experiments reported null results, our analysis implies 
that the either the HI signal is spurious and due to an unknown systematic effect, or 
alternatively that it corresponds to a real resonance, whose production mechanism in DIS 
is substantially different from the production mechanism in e + e~ and the Tevatron. Yet 
another possibility is that either the theoretical or experimental analysis is missing an 
essential ingredient. 

The value of F co ~ 10 required to account for the HI data is surprisingly large. In the 
event that further analysis of the data implies a smaller signal, or an upper limit, the value 
of -F co can be scaled down to obtain the correspondingly reduced predictions for e^e and 
pp. It should be noted that only the overall normalization of the predicted distributions 
and not their shapes will be affected by this procedure. As remarked in sect. 2.2, one 
possible effect that would lead to a reduced estimate of F co is a higher proton yield in 
charm production in DIS, relative to the HERWIG prediction. 
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